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Abstract. This investigation is devoted to the effects of nonadiabatic resonant dynamic tides generated in a 
uniformly rotating stellar component of a close binary. The companion is considered to move in a fixed Keplerian 
orbit, and the effects of the centrifugal force and the Coriolis force are neglected. Semi-analytical solutions for the 
linear, nonadiabatic resonant dynamic tides are derived by means of a two-time variable expansion procedure. 
The solution at the lowest order of approximation consists of the resonantly excited oscillation mode and displays 
a phase shift with respect to the tide-generating potential. Expressions are established for the secular variations 
of the semi-major axis, the orbital eccentricity, and the star’s angular velocity of rotation caused by the phase 
shift. The orders of magnitude of these secular variations are considerably larger than those derived earlier by 


Zahn (1977) for the limiting case of dynamic tides with small frequencies. For a 5 Mq ZAMS star, an orbital 
eccentricity e = 0.5, and orbital periods in the range from 2 to 5 days, numerous resonances of dynamic tides 
with second-degree lower-order (/"''-modes are seen to induce secular variations of the semi-major axis, the orbital 
eccentricity, and the star’s angular velocity of rotation with time scales shorter than the star’s nuclear life time. 


1. Introduction 


Various investigations have already been devoted to tidal 
actions in components of close binaries and their effects 
on orbital evolution. 

A commonly used approach for the study of dy¬ 
namic tides rests on a Fourier decomposition of the tide¬ 
generating potential in terms of multiples of the orbital 
frequency. For each forcing frequency, the dynamic tide 
is then determined by integration of the system of differ¬ 
ential equations governing linear, forced oscillations of an 
equilibrium star on the assumption of a fixed orbit. 


Zahn ( 1975 ) used this approach in order to determine 
the influence of the radiative damping in the nonadia¬ 
batic surface layers of a spherically symmetric star on low- 
frequency dynamic tides. He considered a star consisting 
of a convective core and a radiative envelope. In developing 
an asymptotic theory, he distinguished between the inte¬ 
rior layers, where the forced oscillations are nearly isen- 
tropic, and a boundary layer, where the radiative damping 
is important. 

Because of the radiative dissipation near the star’s sur¬ 
face, dynamic tides do not have the same properties of 
symmetry as the exciting potential, so that the companion 
exerts a torque on the star which may tend to synchronise 
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the star’s rotation with the orbital motion. According to 
Zahn, the torque appears to be strong enough to achieve 
synchronisation, for relatively close binaries, in a time 
short in comparison with the star’s nuclear life time. 
Consequently, Zahn considered dynamic tides with radia¬ 
tive dissipation as an efficient process for the synchroni¬ 
sation of components of close binaries that do not have a 
convective envelope. 


A generalisation of Zahn’s treatment to include the 
effects of the Coriolis forc e in th e case of a slowly rotating 
component was given by Rocca ( 1989 ). 

The system of equations governing the non-adiabatic 
tidal response of a spherically symmetric star was inte¬ 
grated numerically by Savonije & Papaloizou (1983). For 
simplicity, the authors restricted themselves to small or¬ 
bital eccentricities and neglected the perturbation of the 
gravitational potential caused by the tidal response. They 
presented results for main sequence, non-rotating tidally 
distorted stars and concluded that the “the spin up of the 
massive component appears to be rather inefficient ..., 
unless the sys tem could be locked in t o res onance ... ”. 
Subsequently, Savonije fc Papaloizou (|1984| ) determined 
the tidal response for an evolved stellar model of 20 Mq 
and found a substantial increase of the effectiveness of 
tidal evolution during the late phases of core hydrogen 
burning. 

On his side, Tassoul ( 1987 ) argued that the general 
trend toward synchronism that is found in early-type bi- 
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naries, cannot be explained solely by Zahn’s theory on 
the action of radiative damping on dynamic tides. He pre¬ 
sented a hydrodynamical mechanism involving large-scale, 
mechanically driven currents which transport angular mo¬ 
mentum. 


More recently, the work of Savonije & Papaloizou was 
extended by a series of investigations on uniformly ro¬ 
tating spherically symmetric stars in which the effects of 


the Coriolis force were taken into account (Savonije et al. 

1995; 

Papaloizou & Savonije, 

1997; 

Savonije & Papaloizou 

1997; 

Lai, 

1997; 

Witte & Savonije 

1999a). 


The effect of resonance locki ng on the tidal evolution of 
mass ive binaries was studied by Witte fc Savonije ( 1999b , 
2001 ). The authors concluded that when both orbital and 
stellar evolution are taken into account, a dynamic tide 
can easily become locked in a resonance for a prolonged 
period of time. They furthermore found that, during a res¬ 
onance, a significant enhancement of the secular evolution 
of the semi-major axis and the orbital eccentricity occurs. 


In this study, our aim is first to describe the influ¬ 
ence of the perturbation of the radiative energy transfer 
on tidal resonances with free oscillation modes in com¬ 
ponents of close binaries. To this end, we use a two-time 
variable expansion procedure leading to semi-analytical 
solutions for the resonant dynamic tides and the resulting 
perturbation of the star’s external gravitational field. This 
investigation is an extension of an earlier investigation of 
Smeyers et al. (1998) (hereafter referred to as Paper I), 
in which the resonant dynamic tides were treated in the 
isentropic approximation. 


In the second part of our investigation, we determine 
the effects of the nonadiabaticity of the tidal motions on 
the orbital evolution of a close binary in cases of reso¬ 
nances of dynamic tides with free oscillation modes g~^. 


We consider a uniformly rotating component of a close 
binary but we still neglect the effects of the centrifugal 
force and the Coriolis force on the resonant dynamic tides. 

The plan of the paper is as follows. In Sect. 2, we 
present the basic assumptions and the equations governing 
nonadiabatic forced oscillations of a rotating spherically 
symmetric star. In Sect. 3, the tide-generating potential is 
expanded in terms of spherical harmonics and in Fourier 
series in terms of the companion’s mean motion. In Sect. 4, 
we derive the lowest-order solutions for the nonadiabatic 
resonant dynamic tides. Sect. 5 is devoted to the deriva¬ 
tion of the equations governing the secular changes of the 
semi-major axis, the orbital eccentricity, and the longi¬ 
tude of the periastron due to the perturbation of the star’s 
external gravitational field. In Sect. 6, we determine the 
torque exerted by the companion on the tidally distorted 
star and the associated rate of secular change of the star’s 
rotational angular velocity. In Sect. 7, we present phase 
shifts for resonant dynamic tides in a 5 Mq zero-age main 
sequence star. We also consider time scales of orbital evo¬ 
lution for short orbital periods ranging from 2 to 5 days 
and for the orbital eccentricity e = 0.5. Finally, Sect. 8 is 
devoted to concluding remarks. 


2. Basic assumptions and governing equations 


Consider a close binary system of stars with masses Mi 
and M 2 that are orbiting around each other under the in¬ 
fluence of their mutual gravitational force. The first star is 
assumed to rotate uniformly around an axis perpendicular 
to the orbital plane in the sense of the orbital motion, and 
the second star, called the companion, is considered to be 
a point mass describing a fixed Keplerian orbit around the 
star. 

As in Paper I, we describe the oscillations of the 
star with respect to a corotating frame of reference 
whose origin coincides with the mass centre 
Cl of the star and whose a:"^a:"^-plane coincides with the 
orbital plane of the binary. The direction of the a:"^-axis 
corresponds to the direction of the star’s angular velocity 
O. With respect to this corotating frame of reference, we 
introduce a system of spherical coordinates r = 

We consider the spherical coordinates r, 9, (j) as generalised 
coordinates q^. 

The tidal force generated by the companion is intro¬ 
duced as a small time-dependent force which perturbs the 
hydrostatic equilibrium of a spherically symmetric static 
star. Correspondingly, the tide generated by the compan¬ 
ion is considered as a forced perturbation of the star’s 
equilibrium. 

Let StW (r, t) be the tide-generating potential, where 
Et is a small dimensionless parameter defined as 


Here i?i is the radius of the spherically symmetric equilib¬ 
rium star, and a the semi-major axis of the companion’s 
relative orbit. 

Furthermore, let with j = 1,2,3, be the con- 

travariant components of the tidal displacement in the star 
with respect to the local coordinate basis d/dq^, d/dq^, 
d/dq^. If one neglects the effects of the centrifugal force 
and the Coriolis force, these components are governed by 
the equations 
52 {Sq^)j, 

9tj 


+ Utj {5q^)j. = -Et 


dW 


i = 1,2,3. 


( 2 ) 


df^ ' dq^ 

The qij are the components of the metric tensor, and the 
Uij the tensorial operators applying to free linear, isen¬ 
tropic oscillations of a spherically symmetric star, which 
are defined by 




d (<5$).^ (Mt 9P 


1 d {SP)j 
P 


( 3 ) 


dq^ dq^ p dq^ 

In these equations, P is the pressure, p the mass density, 
and $ the gravitational potential. When the operator 5 is 
applied to a quantity, the Lagrangian perturbation of that 
quantity is taken. 

The Lagrangian perturbations of the mass density and 
the gravitational potential are determined by the equation 
expressing the conservation of mass 


= -V, {5q^), 


( 4 ) 
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and by Poisson’s integral equation 

mT = {Sq^)T^j^ 


-G 


[ p{r') {Sq'^)j,{r')[V,j\r -r'\ dV {r'). 
Jv 


(5) 


In the latter equation, G is the gravitational constant, and 
V the volume of the spherically symmetric equilibrium 
star. The operators Vj and are the operators of partial 
differentiation with respect to the generalised coordinates 
and gb ■y^rhen they are applied to a scalar quantity, and 
are the operators of covariant differentiation when they 
are applied to a vector or a tensor component. 

For tides considered as linear, nonadiabatic forced per¬ 
turbations of the spherically symmetric equilibrium star, 
the Lagrangian perturbation of the pressure is determined 
by the equation for the rate of change of thermal energy 


at 


P 


-Ti 


{Sp)j 


P 


(r 3 -i)p 


dt 


( 6 ) 


(7) 


Here Fi = (91nP/cHnp)s and Fa — 1 = {dliiT/dlxip)s 
are generalised isentropic coefficients, T is the tempera¬ 
ture, and dQ/dt is the rate of energy exchange per unit 
time and unit mass. In the regions of the star where the 
energy transfer is purely radiative, dQ/dt is given by 

-TT — -V • Fr, 

dt p 

where £i is the amount of energy generated by the ther¬ 
monuclear reactions per unit time and unit mass, and F r 
the radiative energy flux. 

We make the equations dimensionless by expressing 
the time t, the radial coordinate r, the mass density p, 
the pressure P, both the gravitational potential and 
the tide-generating potential £t W, and the rate of en¬ 
ergy exchange per unit mass dQ/dt respectively in the 

units i?i, Mi/(47ri?f), GMl/{i-KR\), 

GMi/Ri, and Li/Mi, where Li is the luminosity of the 
star. 

Furthermore, we express the Lagrangian perturbation 
of the rate of energy exchange during the tidal motions 
formally in terms of the components of the tidal displace¬ 
ment as 

= (8) 

where the Vj are linear operators. Eq. (^) for the rate of 
change of thermal energy can then be rewritten as 


m 

where 


{SPh 

P 




= c 


(r 3 -i)p 

p 


G = 


GMi 


1/2 


/ 


GM/ 


(9) 


( 10 ) 


Ri Li ^ 

is the ratio of the dynamic time scale to the Helmholtz- 
Kelvin time scale of the star. 


3. The tide-generating potential 

Following Polfliet fc Smeyers ( 1990| ), we expand the tide¬ 


generating potential in terms of unnormalised spherical 


harmonics Y/^{9, (/) and in Fourier series in terms of mul¬ 
tiples of the companion’s mean motion n. The expansion 
takes the form 

A i oo 

Et if (r, t) = Et EE E IF^.m.fe {r ) 

i—2 m——i k— — oo 


exp [i (ctt t — fc n r)], (11) 

where 

We,^,k{r) = -c^,m,k Yri0,c/), ( 12 ) 

ctt is the forcing angular frequency with respect to the 
corotating frame of reference determined as 

CTT = fcn- + TOH, (13) 


and T is a time of periastron passage. The Fourier coeffi¬ 
cients Ci^rn,k are defined as 

e -2 


„ _ |m|)- o|™bn^ 

- //5 , I IM / (^) 


RiV 


1 


1 


- (1 

TT./n 


{£+\m\yG^ (i-e2)^-i/2 

e cosvY~^ cos{k M + mv) dv. (14) 


In this definition, pj"*^(a:) is an associated Legendre poly¬ 
nomial of the first kind, and M and v are respectively 
the mean and the true anomaly of the companion in its 
relative orbit. 

^From Definition ([l^), it follows that the Fourier co¬ 
efficients ce^rn,k depend on the orbital eccentricity. They 
decrease as the multiple fc of the mean motion increases. 
The decrease is slower for higher orbital eccentricities. 
Consequently, the number of terms that has to be taken 
into account in Expansion of the tide-generating po¬ 
tential increases with increasing values of the orbital ec¬ 
centricity. For illustration, the logarithms of the absolute 
values of the second-degree Fourier coefficients C 2 ,m,/c, with 
m = —2,0, 2, are displayed in Fig. as functions of fc for 
the orbital eccentricities e = 0.25 and e = 0.5. 


4. Lowest-order solutions for nonadiabatic 
resonant dynamic tides 


We consider a single partial dynamic tide generated by 
the term 


eTlFf.m.fc(r’) exp[i(crT/-fcnT)] (15) 

of Expansion ( 0 ) of the tide-generating potential, and 
assume that its forcing angular frequency (Tt is close to 
the degenerate eigenfrequency of the star’s free linear, 

isentropic oscillation modes of radial order N that belong 
to the spherical harmonics Y/^{d,(j)) of degree i and the 
various admissible azimuthal numbers m. We define the 
relative frequency difference e as 


£ = - 

<yt,N 

and assume e to be of the order of £t'- 


(16) 


e = /et, /elR. (17) 

In addition, we consider the nonadiabatic effects, 
which are rendered by the right-hand member of Eq. (^ , 
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k 


Fig. 1. Logarithmic representation of the absolute values 
of the Fourier coefhcients C 2 ,m,fc, with m = —2,0,2, as 
functions of k, for the orbital eccentricities e = 0.25 (solid 
circles) and e = 0.5 (open circles). 


to be of the order of Et at least in some region near the 
star’s surface and, for convenience, set 

C = C'eT, C G H. (18) 

We determine the combined effects of the resonance 
of the dynamic tide and the nonadiabaticity of the tidal 
motions by means of an appropriate two-tim e variable ex- 
pansion procedure (Kevorkian fc Cole 1981, Section 3.2.7; 


1991: , Section 4.3.1). In this expansion procedure, the non¬ 


adiabatic terms are evaluated in the quasi-isentropic ap- 
pr oximation. A similar p erturbation procedure was used 
by Terquem et al. ( |l998 ) for the evaluation of the nona¬ 
diabaticity of tidal motions. 

In view of the application of the two-time variable ex¬ 
pansion procedure, we pass on to a new dimensionless time 
variable 

= (19) 

Adopting Et as the small expansion parameter, we intro¬ 
duce a fast and a slow time variable as 

t+ = t*[\^0 (4)] (20) 

and 

i=ETt\ ( 21 ) 

and expand the components of the tidal displacement in 

terms of the star’s linear, isentropic oscillation modes s' 

as 

OO 

s' n—0 

where j = 1, 2,3. The expansions can be restricted to the 
star’s spheroidal modes, since neither the tidal force nor 
the nonadiabatic effects produce any vorticity around the 


normal to the local spherical equipotential surface of the 
equilibrium star. A given spheroidal mode s' is charac¬ 
terised by the degree X' and the azimuthal number fj,' of 
the spherical harmonic to which it belongs, and by the 
radial order ly'. The components of the Lagrangian dis¬ 
placement and the angular eigenfrequency of the oscilla¬ 
tion mode s' obey the wave equations 

4..' {Sq')yyy = 0, * = 1, 2, 3. (23) 

We also introduce an expansion for the Lagrangian per¬ 
turbation of the pressure as 




{r,t+,i) 


(24) 


(SP)^ {r,t) = Y,El 

n—0 

Substitution of Expansions (^) and (^J) into Eq. (^ , 
application of the chain rule, and integration with respect 
to the fast time variable yield, at order e^, 

(5P)® 


P 


= -riE^i' (t+t) 


(25) 

leads 


Introduction of the latter expansion into Eqs. 
to the equations 

s' 

+ ^ 0 (' 5 '?'),. = 0 . * = 1 - 2 , 3 . ( 26 ) 

We transform these equations by using wave Eqs. (^^ and 
multiplying all terms by p{dq‘^)^, where the bar denotes 
the complex conjugate, and s stands for the degree A, 
the azimuthal number and the radial order iz of a given 
spheroidal mode. Next, we integrate over the volume of the 
spherically symmetric equilibrium star and take into ac¬ 
count the orthogonality property of the spheroidal modes. 
The resulting equation for the functions 4°^ (^^>0 takes 
the form 


dt+‘^ 


+ (t+t) = o. 




A general solution is given by 

4°) (4t)=AW (t)expfi^t+) 

V J 


(27) 


(28) 


where (i'j is a yet undetermined complex function of 
the slow time variable t. 

At order et, substitution of Solutions (^ into Eq. (|^ 
and integration with respect to the fast time variable yield 


{6P) 


( 1 ) 


P 


4 °' (t)' 




, (r3-i)p 


x'y 


exp 




4 (V),,. (29) 


In this equality, the second term in the right-hand mem¬ 
ber is the contribution of the nonadiabatic effects to the 
Lagrangian perturbation of the pressure. 

By proceeding in a similar way as for the derivation 
of Eq. (^^, it follows that the functions 4^^ 
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governed by the inhomogeneous equation 


Ul jy 


= -2i- 
1 


dA 


( 0 ) 




dt (Je^N 


, . o-A,, 

exp I 1 ■ 


<yi,N 


a 


l,N 


Afs 


dWe, 


.,k 


s Uv 


dq^ 


m.dV 


I exp [—i (/ f + /c n r) ] | exp (it^) 

iC" 




K 


E 


f - d\(T3-l)pVJSq^) ,] 

[ ^ ^ dV 


4'^ (i) exp [i 

\ J 


(30) 


dq^ 


{6q^)^ dV = -Sx^i5fj.^r. 


2e+l {£-\m\)l 


Ci,m.k / ^ + (^+ l)?7A,j./r]pr^dr, 


(31) 


dAf 

di 


C 


fie,N 


1 


<^e,N 

.,k Qe,N 


a(o) 


2i 


exp [—i (^fi+knr)] . 


(33) 


In this equation, the factor Qi^n is defined as 

1 £r^~^ [^t,7V + (^ + 1) rie,N/r] pr'^ dr 


Qe,N = 


/o”' ^In + + 1) vIn/^'^ 


pr'^ dr 


(34) 


[see Paper I, Eq. (50)], and is the star’s coefficient of 
vibrational stability for the spheroidal mode S defined as 


KiN = — - 


In (^3-l)(dp/p),N [S(dQ/dt)j,^pr^dr 




Jo 


e,N 


+ £(£+ l)p^j^/r^ pr'^dr 


(35) 


(Ledoux & Walraven, 1958, Sect. 63). The star is vibra- 


tionally stable with respect to the spheroidal mode S when 
ki^n is positive, and vibrationally unstable with respect to 
that mode when is negative. 

A general solution of Eq. (^) is given by 


where Afg is the square of the norm of the Lagrangian 
displacement for the oscillation mode s [Paper I, Eq. (48)]. 

The integral contained in the second term of the right- 
hand member of Eq. (^) is developed in Eq. (44) of 
Paper I. After passing on to dimensionless quantities, one 
has 

f dWg^m.k r r ^ TT {£ + \m\)\ 


A 


(0) 


( 0 ) 

Oc exp 


S (^ ) — “s 

— iC l^l,N 


-C 


Ke,N 

<^£,N 


t -b 


£t cke,N 


^ ^e,N 

AO) 


(J2 




ci,m,k Qe,N exp [-i{fi+kn r)] ,(36) 


where <5 a/ is a Kronecker’s delta, and ^a,i/(j') and ri\^^{r) 
are the radial parts of the radial and the transverse compo¬ 
nent of the Lagrangian displacement [Paper I, Eqs. (43)]. 

The integral contained in the third term of the right- 
hand member of Eq. ( p0| ) can be transformed by means of 
Gauss’ theorem. If the mass density vanishes at the star’s 
surface, one derives that 

_ 47 r {X+\p\y. ^ 

2A + 1 {X-\p\)l M.A 

f iT3-l){Sp/p)^^[didQ/dt)]y ^,pr‘^dr. (32) 

^0 

In the inhomogeneous part of Eq. (^0|), the resonant 
terms must be removed. We distinguish between two cases 
depending on whether the star’s spheroidal mode s is 
equal to or different from the spheroidal mode with eigen- 
frequency (Te^N that belongs to the spherical harmonic 
(j)) and is involved in the resonance. We denote the 
latter spheroidal mode as the spheroidal mode S. 

When s = S, the resonant terms are removed from the 
inhomogeneous part of Eq. (^) by setting 


where ag is an undetermined constant. Eor the sake of 
simplification, we set Og ^ = 0. 

Next, when s ^ S, one removes the resonant terms 
from the inhomogeneous part of Eq. ( ^0[ ) by setting 

A® (i) = 0. (37) 

This set of modes s includes the spheroidal modes with 
eigenfrequency ae^N that belong to the spherical harmon¬ 
ics Y^(9, 4>) of degree £ but of azimuthal number p ^ m. 

By introducing the phase angle ipk,s of the complex 
number e — i C K£^n as 
C Kl^pj 


tauA'fe.s = -- 


(38) 


£ Ckl,N 

it follows that the lowest-order solutions for the compo¬ 
nents of the resonant dynamic tide take the form 
£T Ci^m,k Qt,N 


{5q^)^{r,t) = 


1/2 




exp[i(CTTt-fcnr-b'(/’fe.s)], j = 1,2,3. (39) 

Hence, at the lowest order of approximation in the ex¬ 
pansion parameter e^, the solution consists of the reso¬ 
nantly excited oscillation mode with eigenfrequency 
that is associated with the spherical harmonic Yl^{9,(j)). 
This conclusion is similar to the conclusion reached in 
Paper I for resonant dynamic tides considered in the isen- 
tropic approximation. However, the amplitude of the res¬ 
onant dynamic tide is now proportional to 

£t C£,m,kl (e^ + 

and is reduced in comparison to the amplitude reached 
by a resonant dynamic tide considered in the isen- 
tropic approximation. The reduction results from the term 
C'^ N!N denominator and is independent of 

the sign of the star’s coefficient of vibrational stability for 
the oscillation mode S. A reduction of the amplitud e due 
t o the nonadiabatic effects was already observed by Zahn 
(1975) [see the comments below his Eq. (2.40)]. 
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A main difference with the solutions found in the isen- 
tropic approximation is that the nonadiabatic effects in¬ 
duce a phase shift ipk,s between the tidal motions as¬ 
sociated with the resonant dynamic tide and the tide¬ 
generating potential. The resonant dynamic tide lags be¬ 
hind the tide-generating potential when the star is vi- 
brationally stable with respect to the oscillation mode S, 
whereas it precedes the tide-generating potential when the 
star is vibrationally unstable with respect to that mode. 
The existence of a phase shift due to the nonadiabatic- 
ity of the tidal motions was taken into consideration by 
Zahn in his earlier investigation [1975, Eq. (4.5)]. Outside 
resonance, a phase shift does not appear until the next 
order of approximation in the small expansion parameter 


(|Willem4 POOOl ). 


As noted by Terquem et al. (1998) in connection with 
an expression for the torque caused by dissipation in the 
convective envelope of a solar-type star, Solutions (^ are 
valid only for forcing angular frequencies ctt for which 


i^\n/ (40) 

since the nonadiabatic effects are assumed to be small in 
the linear perturbation analysis. 

Besides the partial dynamic tide with forcing frequency 
(Tt, we also consider the partial dynamic tide with forcing 
frequency —ax- The latter partial dynamic tide is resonant 
with the isentropic oscillation mode with eigenfrequency 
One obtains the corresponding solutions for the 
components of the tidal displacement by replacing k by 
—k and m by —m in Solutions (p^. 

When one takes into account the property that 
ci^-m,-k = ci,m,k and reintroduces the dimensions for the 
various physical quantities, the global solutions for the 
components of the resonant dynamic tide take the form 


{5r)T {r,t) = 


P\'^\cose) cos/3m,fe(((),t), 
{56)t {r, t) = -- ^ Qt,N 

m,N{r) (cos 6») 


Et 


(e2 + 


1/2 


t,N 


do 


cos drn^k (0; ^); 


{S(j))T {r,t) = - 


Et Ci,m,k 




1/2 


Qe, 


N 


^^w(e) ^ sin/3m fc(0,t), 

sm 0 


where 


(41) 


f3m,k{4',t) = m(l) +art - kuT + (42) 

Since both the dynamic tide with forcing frequency 
ut and the dynamic tide with forcing frequency —ax are 
taken into account in the derivation of Solutions (^l|), only 
non-negative values of k must be considered subsequently. 


5. Secular variations of the orbital elements 


The tidal distortion of a star brings about a perturba¬ 
tion of the star’s external gravitational field, which in its 
turn gives rise to time-dependent variations of the orbital 
elements. For nonadiabatic resonant dynamic tides, the 
phase shift between the tidal displacement field and the 
tide-generating potential induces variations of the semi¬ 
major axis and the orbital eccentricity besides the varia¬ 
tion of the longitude of the periastron, which is already 
found in the isentropic approximation. 

By means of a classical perturbation procedure used in 
celestial mechanics, the equations governing the rates of 
change of the semi-major axis a, the orbital eccentricity 
e, and the longitude of the periastron w are derived from 
a perturbing function R as 


da 

2 dR 

dt 

v? a dr' 

de 

1 

(1 

dt 

na^ e 


dw 

1 (1- 

-e^ 

dt 

na^ 

e 


n 

\l/2 


dr 

dR 


+ (1 - 


1/2 dR 
dw 


(43) 


de 


(e.g., pternej |l960|) . 

The perturbing function R is related to the Eulerian 
perturbation of the star’s external gravitational potential 
at the instantaneous position of the companion as 


R{u, v,t) = — 


M 1 + M 2 
Wi 




(44) 


where u is the companion’s radial distance. 

In accordance with the lowest-order Solutions (|4^) for 
the components of a nonadiabatic resonant dynamic tide, 
the Eulerian perturbation of the gravitational potential at 
the star’s surface takes the form 

<^T ^ T / / 7-1 \ 

i.Rl) 




(e 2 -h ^IxiI^In) 


Pj’"'(cOs6») cos (3rn,k{4>,t)- 


(45) 


Because of the continuity of the gravitational potential at 
the star’s surface, the Eulerian perturbation of the exter¬ 
nal gravitational potential can be expressed as 


K (r, t) = 


Et Ce^m,k 


Y/J Qi,N ^'i,N (Rl) 


(52 -H 
/ {^+1} 

Pj’”'(cos6»)cos/3m,fc((/),t). 


(46) 


With the use of Definition (|^) of the small parameter ex, 
the perturbing function R can then be written as 

C£ m,k 


R{u^ V, t) = 




1/2 


G (Ml -H M 2 ) p\m.\ 


a 


1+4: 


M 2 


Ml \a 




cos {mv + k M + ij^k.s) ■ (47) 


Here is a dimensionless factor defined by Eq. (62) 
of Paper I. This factor renders the combined effect of the 
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overlap integral, which is proportional to the work done by 
the tidal force, and the Eulerian perturbation of the grav¬ 
itational potential associated with the resonantly excited 
oscillation mode. 

Substitution of the expression for the perturbing func¬ 
tion into Eqs. (^3|), transformation of the time derivatives 
into derivatives with respect to the mean anomaly M, and 
averaging the derivatives da/dM^ de/dM, and dw/dM 
over a revolution of the companion yield the following 
equations for the relative rates of secular change of the 
orbital elements a, e, and vj resulting from a nonadiabatic 
resonant dynamic tide: 

a V Aec V a y Ml Torb 




2 „2 


a: 


e,N 


K 




t,N 


(48) 


1 /de\ _ 27r 

e WVsec V a / Ml Torb 




■Hi, 


G 


(3) 


c(e) 


N 






(49) 


1 / dw\ 

27 ^ V dt 


a 


1.+3 


M. 


Ml Torb 


ecr? 


2^x2 


cr 




(50) 


'i,N ^ ^e,N 

In these equations, Torb is the orbital period, the functions 
^.(e) correspond to the functions Gf_m,fc(e) defined by 
Expression (71) of Paper I, and the functions G^ ^ i^{e) 
and G^^l^ f. (e) are defined as 


G 


( 2 ) 

i,m,k 


(e) = 


:j^Cl^rn,kPi (0)7r 




(£ + 1) e / (1 + e cost?)^ sin(mt? + kM) sinvdv 

Jo 


— m 


/ (1 + e cos cos(m v k M) dv 

Jo 


(51) 


G^iUie) = 


■Q,m.feP’j™'(0)77 


e(l-e2)^ 

|(£-|-l)e J {1 + e cosvY sm{mv + k M) sinv dv 
— mf (1-I-e cosu)^“^ [(1-b e cosw)^ — (l — e^)] 

Jo 


cos{m V + k M) dv j. (52) 

^From Eqs. ( ^ ) - (p0|), it results that the relative rates 
of secular change of the semi-major axis, the orbital ec¬ 
centricity, and the longitude of the periastron due to a 
nonadiabatic resonant dynamic tide are proportional to 
{Ri/aY~^^ and to the mass ratio M 2 /Mi. These factors 
correspond to those found for the rate of secular change 
of the longitude of the periastron in the cases of resonant 
dynamic tides that are treated in the isentropic approxi¬ 


mation. The same factors also appear in various expres¬ 
sions derived for the relative rates of secular change of 
orbital elements outside any resonance with a free oscilla- 


tion mode (Zahc, 

1977, 11978; 

Hut, 

1981; 

Smeyers et al. 

1991; 

Ruymaekers, 

19921). 


Secondly, the relative rates of secular change of the 
semi-major axis, the orbital eccentricity, and the longi¬ 
tude of the periastron depend on the resonantly excited 
oscillation mode through the eigenfrequency the co¬ 
efficient of vibrational stability and the integral ex¬ 
pression Two cases can be distinguished. 

For the relative rates of secular change of the semi¬ 
major axis and the orbital eccentricity, the factor related 
to the resonantly excited oscillation mode is 

—2-;—2— 

^ ^e,N + '^i,N 

The presence of the coefficient of vibrational stability in 
the numerator illustrates that the secular changes of the 
semi-major axis and the orbital eccentricity result from 
the nonadiabaticity of the resonant dynamic tide. 

For the relative rate of secular change of the longitude 
of the periastron, the factor related to the resonantly ex¬ 
cited oscillation mode is 

_ Hi.n/s _ 

1 + kI^/ 

The nonadiabaticity of the oscillation mode thus reduces 
the relative rate of secular change of the longitude of the 
periastron in comparison to the value it takes in the isen¬ 
tropic approximation [Paper I, Eq. (70) 1. T he reduction, 
however, is small because of Inequality (Q). 

Thirdly, the relative rates of secular change of the 
semi-major axis, the eccentricity, and the longitude of the 
periastron depend on the eccentricity through the func- 
tions G(Y^fc(e), Gf)^ ,^{e), and G^Y.fc(e), which are char¬ 
acterised by the degree the azimuthal number m, and 
the number k of the term considered in the tide-generating 
potential. 

For the rates of secular change of the semi-major axis 
and the orbital eccentricity, the order of magnitude in 
the ratio Ri/a can be estimated as follows. When one 
returns to the dimensionless quantities used in the pre¬ 
vious section and expresses the eigenfrequency and 
the coefficient of vibrational stability in the units 

(GMi/i?f)^^^ and RiLi/ (GM^), respectively, the factor 
<^t,N 4 " ^,n) 1 

which appears in Eqs. ( ^^ and (^), takes the form 

G cr£ at j O'! AT + G^ 7c| at) • 

Because of Condition (p0[), the approximation holds 
G (7£,Af 7t£,Af ^ G 
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Since both the relative frequency difference e and the 
nonadiabatic effects are considered to be of the order of 


Et, the factor is of the order of Ej,^ so that 




E'^aj N + 


2 

1,N 


OC 


a 


-3 


(54) 


2 TT 


a 


3/2 


GMi 


Torh 

It follows that 

da\ 


1/2 


M 1 + M 2 
Ml 


1/2 


(55) 


1 


^^y+3/2 


— • ( 56 ) 


T = — 


£T Ci^m,k 


E^ + K 


1,N 




1/2 


(-^ 1 ) ( 


M 2 P, 


m| 






sin (mu + fcM + i/)fe_s). (58) 


T = - 

<Jl,N Ke,N 



/+"M2 , 

/GM2M|\ 

2n 

U J 

' Ml ' 

[ Ml + M 2 ) 

' Torb 


vl/2 


2^2_ 


a: 


■ Hi m G 


l,N 




(4) 

£,m,k 


(e), 


(59) 


d4) 
defined as 


where GT/L j,(e) is a function of the orbital eccentricity 


G^LmAe) =m 


{i + \m\)l 


a 


-{£- 2 ) 


^2 


(60) 


Furthermore, the factor 27r/Torb can be related to the ratio 
Ri/a by means of Kepler’s third law as 


{£ — | to |)! 

The function G^^^ f. (e) is different from zero only for non- 
axisymmetric, nonadiabatic resonant dynamic tides and 
has the same sign as the azimuthal number m. It is re¬ 
lated to the functions g]^}^ ^.(e) and G^fj^ ^.(e) as 

l-e^ 


G^i%,kA = 


(l-e2) 


1/2 


G 


(3) 

£,m,k 


(e)- 


2e 




.(61) 


In particular, for £ = 2, the relative rates of secular change 
of the semi-major axis and the orbital eccentricity are of 
the order of (i?i/a)^^^. This relative rate of secular change 
of the or bital eccentr icity is considerably larger than that 
found by ^ahn| (|l977 ) for the limiting case of small forcing 
frequencies: in that case, the relative rate of secular change 
of the orbital eccentricity is of the order (i?i/a)^^^^. 

Finally, it may be observed that the relative rates of 
secular change of the semi-major axis and the orbital ec¬ 
centricity due to a nonadiabatic resonant dynamic tide 
may vary considerably from one mode to another by the 
value of the factor Hi^m- 


6. Secular change of the star's rotational velocity 

As well because of the phase shift induced by a nona¬ 
diabatic resonant dynamic tide, the companion exerts a 
torque T on the tidally distorted star, which is normal 
to the orbital plane. This torque can be derived from 
Newton’s law of action and reaction as the opposite of 
the torque exerted by the tidally distorted star on the 
companion due to the perturbation of the star’s external 
gravitational field: 

T = M 2 {rx V4>'J (57) 

(tZahnI , |1975|) . 

By the use of Solution (^ ) for the Eulerian perturba¬ 
tion of the external gravitational potential, the torque is 
given by 

■ Qe,N 


The tidal torque T exerted by the companion on the 
tidally distorted star affects the angular velocity of the 
star’s rotation in accordance with the conservation of total 
angular momentum of the binary 

“-¥) ="■ 

where the orbital angular momentum £orb is determined 
as 


^orb 


= / 


gm?m^VG 


, 1/2 


(1 - e^) 


1/2 


(63) 


\Mi+M 2 , 

(see, e.g., Alexander , 19731) . Differentiation with respect 
to time and use of Eqs. (|4^), (|4^), and (^) yield for the 
rate of secular change of rotational angular momentum 


dCr. 


dt 
2 TT 


To, 


■ a 


1/2 , 


A) 

0'1,N Ke,N 
> 9 ^ y 

■ ni 


£+3 


Ml 


fGM? Mi 

VMi 


M 2 


1/2 


■H£,NG<-Ar^^,ie). 


(64) 


t orb ^ ^i,N 

Hence, angular momentum is exchanged between the 
star’s rotation and the orbital motion for nonaxisymmet- 
ric, nonadiabatic resonant dynamic tides. 

Under the assumption that the star rotates as a rigid 
body, the rotational angular momentum is determined as 

Grot=/ri, (65) 

where I is the star’s moment of inertia with respect to the 
rotation axis. Since, at the lowest-order of approximation, 
the rate of secular change of the star’s moment of inertia 
due to a resonant dynamic tide is equal to zero, differen¬ 
tiation of Equality (|^ ) with respect to time and use of 
Eq. (64) yield for the relative rate of secular change of the 


star’s angular velocity 

1 ^dn\ 

~ VV; Ml 

cr£,N Ke,N 


n 


dt 
2 TT 


sec 

1/2 


GM^M^yG 
Ml 
(4) 


M 2 / 


Torb I ^ 


He.NGZJe). 


( 66 ) 


After averaging over a revolution of the companion, one 
obtains 


'^£,N + ^£,N 

One derives the order of magnitude in the ratio Ri/a 
for the relative rate of secular change of the star’s an¬ 
gular velocity by proceeding as for the relative rates of 
secular change of the semi-major axis and the orbital ec¬ 
centricity. By expressing the eigenfrequency the co¬ 
efficient of vibrational stability ki^m, and the star’s mo- 
ment of inertia I respectively in the units (GMi/i?f) , 
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RiLi/ (GMj), and MiR\, setting ^ with 7 G ]R, 

and using Kepler’s third law, it follows that 

Hence, for f = 2, the relative rate of secular change of 
the angular velocity is of the order of which is 

considerably larger than the order of found by 


Zahn (1977) for the limiting case of small forcing frequen¬ 


cies. 


With regard to synchronisations of components of close 
binaries, an appropriate quantity is the relative rate of 
secular change of the difference between the star’s angular 
velocity and the companion’s mean motion: 


1 


fl- 


d{fl — n) 


dt 


Ri 


1+3 


M 2 2 TT 
Ml Torb 




rr 2 


■ He,N 


1 




GM 1 M 2 (4) , , 

nah 


Q — n 
3 


1 ( 2 ) 


2 


(e) 


( 68 ) 


7. Nonadiabatic tidal resonances in a 5Mq 
zero-age main sequence star 

We have determined phase shifts for nonadiabatic 
dynamic tides in resonance with lower-order second-degree 
^'''-modes in a 5 Mq zero-age main sequence stellar model. 
The model consists of a convective core and a radiative 
envelope and has a central hydrogen abundance = 0.7 
and a radius i?i = 2.8 Rq. The dynamic time scale of the 
model is 54.8 minutes, the Helmholtz-Kelvin time scale 
4.88X lO® years, and the nuclear time scale 8.65 x lO^ years. 
The ratio C of the dynamic time scale to the Helmholtz- 
Kelvin time scale is equal to 2.14 x 10“^°. 

We have considered resonances with the second-degree 
g+-modes of radial orders N = I, 2,..., 20. The eirenfre- 
quencies ( 72 ,Af of these modes are given in Table ^ For 
the determination of the coefficient of vibrational stability 
within the framework of the quasi-isentropic approxima¬ 


tion, we adopted the expressions given by Degryse et al. 


(1992). Since the quasi-isentropic approximation breaks 
down near the star’s surface, the integration in the nu¬ 
merator of Definition ( ^5| ) for the coefficient of vibrational 
stability is often stopped at the transition between the 
quasi-isentropic stellar interior and the nonadiabatic sur¬ 
face layers (see, e.g., Stothers fc Simon, 197C; Cox, 19801). 


We therefore determined the numerator in Definition 
by breaking off the integration at various radial distances 
close to the star’s surface and found the value of the inte¬ 
gral to be rather insensitive to the radial distance at which 
the integration was stopped. However, the resulting val¬ 
ues of the coefficient of vibrational stability K 2 ,n differed 
strongfy from, or were even in contradiction with results of 
common nonadiabatic stability analyses of slowly pulsat¬ 


ing B stars ( see, e.g., [Dziembowski et al.| , |1993| ; [Gautschy 


fc Saioj |1993| ). Several modes that were expected to be un¬ 
stable turned out to be stable within the framework of the 


Table 1. The eigenfrequencies and the linear growth rates 
for the second-degree modes - g^Q of the 5 Mq ZAMS 
model. 


mode 

<72,iV (s“^) 

«2,iV (S“^) 

dt 

0.661 10“® 

0.372 10"'^ 

at 

0.460 lO-’^ 

0.127 10"'^ 

at 

0.350 lO”’’ 

0.552 10"^ 

at 

0.281 10“® 

0.533 10"'^ 

at 

0.234 10“® 

0.325 10"’^° 

at 

0.200 10“® 

0.315 10"“ 

at 

0.175 10“® 

0.168 10"’^° 

at 

0.156 10“® 

0.182 10"’^° 

at 

0.141 10“® 

0.224 10"“ 

ato 

0.129 10“® 

-0.350 10"“ 

ati 

0.119 lO"’^ 

-0.126 10“® 

ati 

0.110 10“® 

-0.220 10"® 

ata 

0.102 lO"’^ 

-0.271 10“® 

ati 

0.955 10“^ 

-0.117 10“® 

ats 

0.897 10“^ 

0.489 10"'^ 

ate 

0.846 10-* 

0.175 10"® 

atr 

0.800 10-* 

0.391 10"® 

ats 

0.758 10“^ 

0.741 10"® 

ata 

0.720 10“^ 

0.127 10"^ 

ato 

0.686 10“'‘ 

0.201 10"^ 


quasi-isentropic approximation, and vice versa. Therefore, 
we have not used the coefficients of vibrational stability 
determined by means of the quasi-isentropic approxima¬ 
tion but adopted the linear growth rates, which are the 
opposites of the imaginary parts of the complex eigenfre¬ 
quencies of the nonadiabatic modes depending on time by 
exp[i(( 7 -f i«;)f]. This change is justified by the fact that a 
first-order perturbation analysis of resonant dynamic tides 
within the framework of a fully nonadiabatic treatment 
yields solutions for these tides that are formally identical 
to Solutions (^^, except that the linear growth rate of the 
mode appears instead of the coefficient of vibrationaf sta¬ 
bility, and that the factor is determined in terms 

of the nonadiabatic eigenfunctions (see the Appendix). 
Since, for main-sequence stars, the real parts of the nona¬ 
diabatic eigenfunctions are known to be much larger than 
the imaginary parts and almost equal to the isentropic 
eigenfunctions (see, e.g., Unno et al., 1989), the nonadia¬ 


batic factor Qe^N has been approximated by means of the 
expression given by Equality (Q). 

The linear growth rates K 2 ,n for the 5 Mq ZAMS stel¬ 
lar model are presented in Table [^. The model is vibra- 
tionally stable for all 5 +-modes considered, except for the 
modes g^Q - 5 + 4 . 


We have determined the phase shifts on the sup¬ 
position that e = 0.001. The resulting values are displayed 
in Fig. H as a function of the radial order of the mode. 
Vibrationally stable modes are represented by solid cir¬ 
cles, and vibrationally unstable modes by open circles. 

In addition to the phase shifts, we have determined 
the characteristic time scales ta, G, and Gyn associ¬ 
ated with the relative rates of secular change of the semi- 
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Fig. 2. Variation of the logarithm of the phase shift 
'4’k,s as a function of the radial order N of the mode 
for the 5 Mq ZAMS stellar model, and for e = 0.001. 
Vibrationally stable modes are represented by solid cir¬ 
cles, and vibrationally unstable modes by open circles. 



major axis, the orbital eccentricity, the angular velocity, 
and the difference between the star’s angular velocity and 
the companion’s mean motion as 


1 

ta 

1 

tn 


I /da\ 1 

« WVsec ’ te 


1 

l'dQ\ 

n 



1 /de\ 




e \dt) 

sec 

5 


I 

d(f 2 — n) 


— n 


dt 

sec 


For this, we have set M 2 = 1.4 Mq and have adopted a 
value for the star’s rotational angular velocity equal to 
50% of the companion’s orbital angular velocity at the 
periastron. Furthermore, we have considered the relatively 
large orbital eccentricity e = 0.5. 

In the determination of the time scales, we have taken 
into consideration a large number of terms in the expan¬ 
sion of the tide-generating potential. We have also taken 
account of the limitations imposed by the perturbation 
theory: since e is assumed to be of the order of Et and 
must satisfy Inequality we have restricted the values 
of e to the range of values extending from e = O.Iet to 
e = 10 Et- 

The contributions of the nonadiabatic resonant dy¬ 
namic tides to the secular variations of the orbital ele¬ 
ments and of the star’s angular velocity have been re¬ 
stricted to the dominant contributions, which are those 
associated with the azimuthal number m = —2. 

For orbital periods ranging from 2 to 5 days, the vari¬ 
ations of the time scales, expressed in years, are displayed 
in Fig. ||on a logarithmic scale. The gaps in the curves ap¬ 
pear at orbital periods at which the relative frequency dif¬ 
ference e is outside the admitted range of values. The black 


lines correspond to orbital periods at which da/dt < 0 , 
de/dt < 0, dVL/dt > 0, and [1/(11 — n)][d(n — n)/dt] < 0, 
and the grey lines to orbital periods at which da/dt > 0, 
de/dt > 0 , dfl/dt < 0 , and [ 1/(11 — n)][d(fl — n)/dt] > 0 . 
The dotted horizontal line represents the logarithm of the 
star’s nuclear time scale. 

The most striking result is that numerous resonances 
occur which have strong effects on the relative rates of 
secular change of the orbital elements and the star’s ro¬ 
tational angular velocity. For the shorter orbital periods, 
many resonances lead to characteristic time scales that are 
considerably shorter than the star’s nuclear time scale. 

The resonantly excited oscillation modes in the con¬ 
sidered range of orbital periods are the g^-modes of ra¬ 
dial orders N = 1,...,12. At several orbital periods, 
even resonances with two or more oscillation modes oc¬ 
cur. The feature near the orbital period of 3.25 days, for 
instance, is caused by resonances with the modes , 

and g^. Furthermore, resonances with vibrationally un¬ 
stable modes contribute to secular changes of the orbital 
elements and of the star’s angular velocity in the sense op¬ 
posite to that of the contributions of the resonances with 
vibrationally stable modes, so that they counteract the 
effects of the latter resonances. At the orbital period of 
3.66 days, for instance, the resonance with the mode g/j, 
contributes to secular changes of the orbital elements and 
of the star’s angular velocity in the sense opposite to that 
of the contributions stemming from the resonances with 
the modes g^ and gf. 

With the exception of the shortest orbital periods con¬ 
sidered, the time scales for the relative rates of secular 
change of the orbital elements and the star’s rotational 
angular velocity are mostly longer than the star’s nuclear 
time scale, in accordance with the results o f earlier inves¬ 
tigations (e.g., Savonije fc Papaloizou , 1983| ). The shortest 
time scales found for the relative rates of secular change 
are about five orders of magnitude smaller than the nu¬ 
clear time scale. With regard to the synchronisation of 
the star’s rotation with the orbital motion of the compan¬ 
ion, the time scales found are generally shorter than the 
time scales associated with the relative rates of secular 
change of the orbital eccentricity. Synchronisation of the 
star’s rotation with the orbital motion of the companion 
at the periastron of its relative orbit can therefore be ex¬ 
pected to occur before the circulari sation of the orbit. A 
similar conclusion was reached by Zahc (1977) for low- 
frequency dynamic tides, after averaging over the effects 
of resonances. 

When the orbital period becomes larger, nonadiabatic 
resonant dynamic tides have a decreasing influence on the 
time scales for the relative rates of secular change of the or¬ 
bital elements and the star’s angular velocity. The under¬ 
lying reason is the rapid decrease of the quantity 772 ,jv for 
higher-order g+-modes, which are the relevant modes for 
resonances of dynamic tides in close binaries with longer 
orbital periods (see also Paper I). 

On the basis of these results, we suggest that nonadi¬ 
abatic resonant dynamic tides may be an efficient mech- 
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Fig. 3. The logarithms of the time scales ta, tg, tn, and tsyn as functions of the orbital period for the 5 Mq ZAMS 
model and for an orbital eccentricity e = 0.5. The black lines correspond to orbital periods at which da/dt < 0, 
de/dt < 0, dVL/dt > 0, and — n)][(i(r2 — n)/dt\ < 0, and the grey lines to orbital periods at which da/dt > 0, 

de/dt > 0, dfl/dt < 0, and [l/(n — n)][d(n — n)/dt] > 0. The dotted horizontal line represents the logarithm of the 
star’s nuclear time scale. 


anism for the circularisation of orbits in close X-ray bi¬ 
naries, especially when the orbital period is short. X-ray 
binaries are thought to have initially highly eccentric or¬ 
bits as a result of the supernova explosion which generates 
the compact star. However, the orbits of several X-ray bi¬ 
naries with short orbital periods have been reported to 
be circular. Two examples are the X-ray binaries Cen X-3 
and SMC X-1. The orbital period of the system Cen X- 
3 is 2.1 days, and that of the sys tem SMC X-1 3.9 days 
( Verbunt fc Van den Heuvel , 1995 ). For close binaries with 
such short orbital periods and highly eccentric orbits, res¬ 
onances of dynamic tides with low-order g^-modes are 
quite possible. 


8. Concluding remarks 

We have studied the effects of nonadiabatic resonant dy¬ 
namic tides in a component of a close binary on the orbital 
elements and on the star’s axial rotation. The compan¬ 
ion is considered to move in a fixed Keplerian orbit, and 
the star is supposed to rotate uniformly around an axis 
perpendicular to the orbital plane, but the effects of the 
centrifugal force and the Coriolis force are neglected. The 
relative difference between the frequency imposed by the 
companion and the star’s eigenfrequency is assumed to be 
of the order of the ratio of the tidal force to the gravity at 
the star’s equator. Furthermore, the nonadiabatic effects 
are considered to be of the same order of magnitude at 
least in some region near the star’s surface. 

We have derived semi-analytical solutions for the lin¬ 
ear, nonadiabatic resonant dynamic tides by applying a 
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two-time variable expansion procedure, in which the ratio 
of the tidal force to the gravity at the star’s equator is 
adopted as the small expansion parameter. In this pertur¬ 
bation method, the nonadiabatic effects are evaluated in 
the quasi-isentropic approximation. At the lowest order of 
approximation, the resonant dynamic tide corresponds to 
the free oscillation mode involved in the resonance. This 


conclusion is similar to the conclusion reached by Smeyers 


et al. ( .998) for resonant dynamic tides treated in the isen- 
tropic approximation. 

Due to the nonadiabatic effects, the amplitude of the 
resonantly excited oscillation mode is reduced in compar¬ 
ison to the amplitude it reaches in the isentropic approxi¬ 
mation. The magnitude of the reduction is determined by 
the ratio of the coefficient of vibrational stability to the 
eigenfrequency of the oscillation mode. 

The main difference with the dominant solutions de¬ 
rived in the isentropic approximation is that the nonadi¬ 
abatic effects induce a phase shift between the tidal dis¬ 
placement and the tide-generating potential. The magni¬ 
tude and the sign of the phase shift depend on the co¬ 
efficient of vibrational stability and on the relative dif¬ 
ference between the forcing angular frequency and the 
eigenfrtjqueiicy of the oscillation mode involved in the res¬ 
onance. The tidal displacement field lags behind the tide- 
generating potential when the star is vibrationally stable 
with respect to the oscillation mode, whereas it precedes 
the tide-generating potential when the star is vibrationally 
unstable with respect to the oscillation mode. 

The phase shift leads to variations of the semi-major 
axis and the eccentricity of the companion’s relative or¬ 
bit, in addition to the variation of the longitude of the 
periastron already found in the isentropic approximation 
of a resonant dynamic tide. The contributions to the rela¬ 
tive rates of secular change of the semi-major axis and the 
orbital eccentricity are determined by Eqs. (^ ) and 
and are of the order of which is much larger 

than the order of derived by Zahn (1977) for 

the limiting case of dynamic tides with small frequencies. 

As well because of the phase shift, the companion ex¬ 
erts a torque on the tidally distorted star, which is nor¬ 
mal to the orbital plane. The average value of the torque 
is given by Eq. For nonaxisymmetric resonant dy¬ 

namic tides, angular momentum is exchanged between the 
orbital motion and the star’s rotation. On the assumption 
that the star rotates as a rigid body, the relative rate of 
secular change of the angular velocity due to a resonant 
dynamic tide is determined by Eq. (|^ and is of the order 
of which is also considerably larger than the or¬ 
der of derived by Zahn (1977) for the limiting 

case of dynamic tides with small frequencies. 

We have applied our results to resonances of dynamic 
tides with second-degree, lower-order g^-modes in a 5 Mq 
zero-age main sequence stellar model consisting of a con¬ 
vective core and a radiative envelope. For the orbital ec¬ 
centricity e = 0.5, numerous resonances of dynamic tides 
with g'''-modes may occur in the range of orbital periods 


from 2 to 5 days. At several orbital periods, even reso¬ 
nances with two or more (;''"-modes may take place. 

Since the quasi-isentropic approximation led to inade¬ 
quate results for the coefficient of vibrational stability, we 
passed on to the linear growth rates determined within 
the framework of a fully nonadiabatic analysis. 

Relative to the resonances, we have determined time 
scales associated with the rates of secular change of the 
semi-major axis, the orbital eccentricity, and the star’s 
rotational angular velocity. The shortest time scales are 
found for the shortest orbital periods and for resonances 
with the lowest-order g'^-modes. They are about five or¬ 
ders of magnitude smaller than the star’s nuclear time 
scale. The time scales for synchronisation of the star’s ro¬ 
tation with the companion’s orbital motion at the peri¬ 
astron in the relative orbit are generally shorter than the 
time scales for circularisation of the orbit. 

Because of the shortness of the time scales, the sys¬ 
tem may be expected to rapidly evolve away from any 
resonance between a dynamic tide and a free oscillation 
mode, so that only a limited amount of angular momen¬ 
tum can be exchanged between the star’s rotation and 
the orbital moti o n of the companion. However, Witte fc 
Savonij^ (19991:, 2001) have shown that when both stel¬ 


lar and orbital evolution are taken into account, a dy¬ 
namic tide can easily become locked in a resonance for a 
prolonged period of time. In these circumstances, the ex¬ 
change of angular momentum might be sustained through¬ 
out the locking of the resonance. 

For conclusion, in components of close binaries with 
short orbital periods and highly eccentric orbits, various 
resonant dynamic tides can be excited which, through 
their nonadiabatic character, may act as an efficient mech¬ 
anism for the synchronisation of the rotation and the cir¬ 
cularisation of the orbital motion of stellar components in 
time scales shorter than their nuclear time scales. 
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Appendix A: Resonant dynamic tides within a 
fully nonadiabatic framework 

In this appendix, we briefly describe how the lowest-order 
solution for a nonadiabatic resonant dynamic tide can be 
derived within a fully nonadiabatic framework. 

Following Van Hoolst ( |l994 ), we express the equations 
that govern linear, nonadiabatic stellar oscillations in the 
form of a first-order, linear differential equation in the vec¬ 
tor z — d^/dt, 6S), where ^ is the Lagrangian displace¬ 
ment, and 6S the Lagrangian perturbation of the entropy. 
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Adding the tidal force exerted by a companion star, one 
has 
dz 


dt 


= Mz + AT, 


(A.l) 


where N = (0, —0), and M is the integro- 
differential operator governing free nonadiabatic stellar 
oscillations as defined in Van Hoolst ( |l994 ). The vector 
z is expanded in terms of linear, nonadiabatic eigenvec¬ 
tors e„ as 


z(t,r) = y^a„(t)e„(r). 


(A.2) 


where the eigenvectors e„ depend on time by exp(cr„t) and 
obey the equation 

Me„ = cr„e„. (A.3) 

The complex angular frequency tT„ can be expressed as 

0'n=Kn+if^n, (A.4) 

where is the oscillation frequency, and the opposite 
of the linear growth rate. 

In order to derive a first-order expression for reso¬ 
nant dynamic tides, we use a multiple time scales method 


(Nayfeh, 1973, 1981). Therefore, we introduce time vari¬ 
ables To and Ti as 

To =t, Ti = ert, (A.5) 

and expand the amplitude as{t) of the resonant mode in 
terms of the small parameter et as 


( 0 ) 

as = 0-8 


ExOg . 


(A.6) 


In addition, we assume the linear growth rate to be small 
and of the order of Et- In contrast with the procedure 
followed in the quasi-isentropic treatment, we here do not 
consider the nonadiabatic effects to be small everywhere 
in the star. 

The application of the method of the multiple time 


scales as described in Van Hoolst (1994) leads to an equa¬ 
tion for the first-order real part of the eigenfunction that 
has the same form as Eq. (^), except for the fact that 
now stands for the linear growth rate, and the fact 
that the factor Qe^N is defined in terms of the complex 
nonadiabatic eigenfunctions of the mode. 
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